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ABSTRACT 

We present the comparison of the three most important ice constituents (water, CO and CO2) 
in the envelopes of massive Young Stellar Objects (YSOs), in environments of different metal- 
licities: the Galaxy, the Large Magellanic Cloud (LMC) and, for the first time, the Small 
Magellanic Cloud (SMC). We present observations of water, CO and CO2 ice in 4 SMC and 
3 LMC YSOs (obtained with Spitzer-IRS and VLT/ISAAC). While water and C0 2 ice are de- 
tected in all Magellanic YSOs, CO ice is not detected in the SMC objects. Both CO and C0 2 
ice abundances are enhanced in the LMC when compared to high-luminosity Galactic YSOs. 
Based on the fact that both species appear to be enhanced in a consistent way, this effect is 
unlikely to be the result of enhanced CO2 production in hotter YSO envelopes as previously 
thought. Instead we propose that this results from a reduced water column density in the en- 
velopes of LMC YSOs, a direct consequence of both the stronger UV radiation field and the 
reduced dust-to-gas ratio at lower metallicity. In the SMC the environmental conditions are 
harsher, and we observe a reduction in CO2 column density. Furthermore, the low gas-phase 
CO density and higher dust temperature in YSO envelopes in the SMC seem to inhibit CO 
freeze-out. The scenario we propose can be tested with further observations. 

Key words: astrochemistry - circumstellar matter - galaxies: individual (LMC, SMC) - 
Magellanic Clouds - stars: formation - stars: protostars. 



1 INTRODUCTION 

The most abundant astronomical ices are water, CO and CO2. 
These ices involve abundant metals that, as gas-phase atoms 
and molecules, are prime coolants of the molecular clouds (e.g. 
lOmukai et al.l l20ich . Efficient cooling during the early collapse 
stages is paramount to allowing star formation to proceed. By 
changing the chemical balance in the molecular cloud, interstel- 
lar ice chemistry is thus an important ingredient in the formation of 
stars from contracting cores. 

In spite of its ubiquity and significance as a cornerstone of in- 
terstellar chemistry, the formation of water is poorly understood. In 
comparison to gas-phase water abundance, the observed water ice 
abundance is too high to have formed by direct accretion from the 
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gas-phase only ^Roberts fe^ Herbst 2002). Therefore surface reac- 
tions on cold dust grains likely play a significan t role in the forma- 
tion of water molecules ( Tielens & Hagen 1982). 

Water ice seems to form very easily in molecular clouds, as it 
requires relati vely low gas densities a nd only moderately low tem- 
peratures (e.g. lHollenbach et al.l2009l) . CO, the most abundant gas- 
phase molecule besides molecular hydrogen, forms exclusively in 
the gas phase. It is thought that a layer of water ice (with impuri- 
ties of CO and other molecules) forms on cold dust grains. During 
pre-stellar collapse the increased density causes virtually all gas 
species to accrete o nto the icy grains, creating a layer dominated 
mainly by CO (e.g. lAikawa et alj|200lh . Thus the CO ice profiles 
towards quiescent regions are dominated by pure CO or CO-rich 
(apolar) matrices, with an increasing contribution of C O in water- 
rich ( or polar) matrices towards YSO environments (e. g.lGibb et al. 
2004). CO2 forms exclusively via surface reactions jMillar et al.l 
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Figure 1. Optical depth spectra of the Magellanic YSOs. The standard star 
spectrum (dotted line) shows the telluric line at ~ 4.68 /im that affects the 
CO ice feature. Dashed lines show the column density integration intervals. 
No CO ice is detected towards the SMC. The water spectra of the SMC 
objects have been smoothed in this figure. 

Il99lh and the observed profiles suggest that these reactions occur 
both in the CO- and water-dominate d layers, the latt er being the 
dominant component (> 50%, Gera kines et al.|[l999h towards all 
environments. 

In this letter we investigate the effect of metallicity on the 
abundance of the three most important contributors to interstel- 
lar ices: water, CO and CO2. Our research is based on observa- 
tions in the Magellanic Clouds (MCs) which ha ve subsolar metal- 
licity (Zlmc ~ \Zq and Zsmc ~ e.g. lRussell & BesseT3 
1989). We present the first observations of CO and CO2 ice to- 
wards Young Stellar Objects (YSOs) in the SMC. We compare the 
ice abundances towards luminous YSOs in the Galaxy, LMC and 
SMC, and relate these to the observed properties of the Interstellar 
Medium (ISM), e.g. gas-phase abundances, dust temperature and 
radiation field strength. We propose a simple schematic scenario 
that can explain the observed relative ice column densities and the 
non-detection of CO ice in the SMC environments. 



2 SPECTROSCOPIC OBSERVATIONS 

The sample of 7 M agellanic YSOs (Tab lep} includes 3 LMC 
object s described by I Oliveira et aD [2009) and Shimonishi et al. 
feOlCh and 4 SMC s ources. IRAS 01042-7 215 in the SMC was 
identified as a YSO bv lvan Loon et al.l d2008h . The three other SMC 
sources are spectroscopically identified as YSOs for the first time 
in this work, and are amongst the brightest of a sample of YSO 
candidates in the SMC (Oliveira et al. in preparation), identified at 



70 fj,m in the S 3 MC survey dBolatto et al]|2007l) . The spectral en- 
ergy distributions of these YSOs suggest that they are young and 
embedded, and of similar total luminosity (a few 10 4 L©), In terms 
of luminosity, the MC sample is comp arable to the Gala ctic sample 
of massive YSOs observed with ISO dGibb et alJliooi. and refer- 
ences therein). 

M-band spectra were obtained with the Infrared Spectrom- 
eter And Array Camera (ISAAC) at the European Southern Ob- 
servatory (ESO)'s Very Large Telescope (VLT) on the nights of 
4 and 5 November 2009 (ESO Programme 084.C-0955, P.I. J.M. 
Oliveira), to search for the CO ice feature at A = 4.67 /im. The 
low-resolution mode of ISAAC was used with a l"-wide slit, re- 
sulting in a spectral resolution of AA ~ 0.012 /j,m around the CO 
feature, sampled by 0.0014 /im pixels. Exposure times varied be- 
tween 45 and 90 min. The standard IR technique of chopping and 
nodding was employed to cancel the high background. Spectra of 
standard stars were obtained and were used to cancel telluric ab- 
sorption in the YSO spectra; telluric lines were used to calibrate 
the wavelength scale to an accuracy AA ~ 0.003 /im. 

The L-band spectra of three of the SMC sources were ob- 
tained with ISAAC at the ESO/VLT on the nights of 28 and 29 Oc- 
tober 2006 (ESO Programme 078.C-0338, P.I. J. M. Oliveira), the 
spectr um of IRAS 01042— 7215 was published by Ivan Loon et al.l 
(2008). The L-band spectra were obtained and reduced in the same 
way as the M-band spectra. The resolving power is A/AA ~ 500. 
Exposure times varied between 60 — 105 min. The hydrogen lines 
in the standard stars left remnants of at most a few per cent of the 
continuum level. Telluric lines were used to calibrate the wave- 
length scale to an accuracy AA ~ 0.002 /im. The 3—4 /im spectra 
of the L MC sources were obtain ed with the AKARI satellite (pub- 
lished bv lShimonishi et alfeOld) . 

All sources were observed with the InfraRed Spectrograph 
(IRS) on board the Spitzer Space Telescope, either as part of the 
SAGE -Spec Spitzer Legacy P rogram of the LMC l lKemper et all 
l201Ch and already described in Oliveira et al. (2009), or as part of 
GTO programme #50240 on the SMC (P.I. G.C. Sloan) and shown 
here for the first t ime. Detai l s on th e data reduction of the IRS spec- 
tra are found in lSloan et all d2008l) . 

Continuum subtraction was performed by fitting a 2 nd - to 5 th - 
degree polynomial over a narrow wavelength range surrounding the 
ice features. Fig. Q] shows the optical depth spect ra. We detect wa- 
ter, CO and CO2 ices towards all 3 LMC sources; Shimonishi et al. 
feOlCh did not detect CO ice towards IRAS 05240-6809. We de- 
tect water and CO2 towards the SMC sources (only at a 3-cr level 
for IRAS 01042— 7215); the most striking observation is the ab- 
sence of CO ice in any of the SMC spectra (Fig.Qjmiddle). 



3 RESULTS: ICE COLUMN DENSITIES 

We measured column densities for the three ice species by in- 
tegrating the optical depth spectra over fixed wavelength ranges 
(2.9-3.6, 4.66-4.69 and 14.8- 16.2 /im). Fig.[[] shows that we 
miss the extreme red wing of the H2O ice profile, which accounts 
for ~10% of the total column density. The integration interval for 
the CO ice excludes the Pf/3 hydrogen recombination line at A = 
4.65 fim and the S(9) line from molecular hydrogen at A = 4.687 
/im. Adopted band strengths are 2.0 x 10~ 16 , 1.1 x 10~ 17 and 
1.1 x 10~ 17 cm per m olecule, respectively for H2O, CO and CO2 
dGerakines et at J 1995 ). For the LMC objects, we compare our mea- 
surements to those by Shimonish i et alj (1 20 lOh . Our water and CO2 
column densities are generally smaller (by 10—20%), due to dif- 
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Table 1. Target positions and ice column densities (in 10 17 molecule cm 2 ). 



Galaxy 


Target ID 




RA & Dec (J2000) 


N(H a O) 


N(CO) 


N(C0 2 ) 


LMC 


SSTISAGE1CJ051449.41 


-671221.5 


05:14:49.41 


-67:12:21.5 


19.3 ± 1.0 


1.60 ±0.11 


6.10 ±0.45 


LMC 


IRAS 05240-6809 




05:23:51.14 


-68:07:12.4 


27.5 ± 1.5 


2.62 ± 0.22 


12.92 ± 0.26 


LMC 


SSTISAGE1C J052546.49 


-661411.3 


05:25:46.52 


-66:14:11.3 


33.6 ± 1.4 


1.54 ±0.16 


8.93 ±0.18 


SMC 


IRAS 00430-7326 




00:44:56.30 


-73:10:11.8 


17.7 ±0.7 


^0.26 


2.08 ± 0.09 


SMC 


S3MC 00540-7321 




00:54:02.31 


-73:21:18.6 


21.6 ±0.8 


^0.19 


3.50 ±0.19 


SMC 


S3MC00541-7319 




00:54:03.36 


-73:19:38.4 


22.3 ± 1.2 


^0.16 


5.40 ±0.16 


SMC 


IRAS 01042-7215 




01:05:49.29 


-71:59:48.8 


16.6 ± 0.7 


^0.45 


0.49 ±0.18 
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Figure 2. Ice column densities. Galactic data (triangles) are from iGerakines etal] ll999l).lGibb et al. (2004) and Pontoppidan et al. 12008, low-luminosity 
YSOs, smaller triangles). LMC data (circles) are from this work and lShimonishi et alj |201CI) ; SMC data (diamonds) are from this work. The solid lines in 
panels a, b and c are line fits in linear space. The ternary plot for the column densities is also shown (panel d, see text); in this diagram the lines associated 
with the CO non-detections (open symbols) indicate the change in object position if the CO column density is a third of its upper limit. 



ferent continuum and integration interval choices. For CO ice, our 
column densities are lower by as much as 50%. This is likely due 
to the lower resolution of the AKARI spectra (R ~ 80, compared to 
R ~ 390 for our spectra) that makes continuum determination un- 
certain. Such effects are not accounted for in the statistical errors. 

Fig-H] shows t he ice column densitie s for the MC samples 
(from this work and Shimonishi et al. 2010) and the massi ve YSO 
Galactic samples dGerakines et alJI 19991 ; fGibb et al.ll20ol) for the 
3 ice species analysed here — we only consider objects with mea- 
surements or upper limits for all 3 species. The N(C02)/N(H20) 
ratios for the LMC and the Galaxy are consistent with previ- 
ous determinations (Gerakines et al. 1999]; loiiveira et all 120091 : 



IShimonishi et alfeOld) . respectively ~ 0.32 and ~ 0.2 (Fig.|2](a)). 
Thus CO2 ice seems to be more abundant (with respect to water 
ice) in the LMC than in the Galaxy. For the SMC sources, this 
ratio is lower than those of the other samples: in comparison to 
both LMC and Galactic objects with similar water column den- 
sity (16.6 -22.3 x 10 17 molecules cm" 2 ), the SMC sources exhibit 



smaller CO2 ice column density. The sample size and water col- 
umn density range sampled are however small. Nevertheless, this 
suggests that while CO2 ice abundance is enhanced in the LMC, it 
is depleted in the SMC. 

CO ice is detected for most of the LMC sources. In Fig.|2] 
(c) CO and CO2 column densities are plotted against each other. 
The LMC and Galactic sources are indistinguishable in this dia- 
gram, and there is a strong correlation (Pearson coefficient ~ 0.94). 
With the exception of a single source, no CO ice is detected for 
sources with CO2 column density < 5xlO lr molecule cm -2 . 
This is expected since CO is the most volatile of the ice species 
considered here. The slope of the linear relationship is ~ 0.84, 
and a CO2 column density threshold for the detection of CO of 
~ 4.7 xlO 17 molecule cm " 2 is inferred (similar to those derived 
bv IGerakines e 7aT] [T999T) . Since for a given CO2 ice abundance, 
Galactic and LMC objects show similar CO ice abundances, this 
means that both species are affected in a similar way by whatever 
causes the CO2 enhancement with respect to water ice. 
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CO and water column densities are clearly correlated for the 
LMC and Galactic samples (Fig.[2](b), Pearson correlation coeffi- 
cient of ~0.8), but there is considerable scatter. Even though less 
clear than in the case of CO2, LMC sources tend to exhibit larger 
CO column densities than Galactic sources for the same water col- 
umn density. Due to the larger scatter, we do not fit these data; 
instead we derive the linear relations that are consistent with those 
observed between CO and CO2, and CO2 and water. This results in 
slopes of 0.27 and 0.17 respectively for the LMC and the Galaxy, 
and water column density thresholds for the presence of CO of ~ 
14.5 and 23 x 10 17 molecule cm" 2 . 

Since water, CO an d CO2 are the most abundant ice species 
(e.g. Ivan Di shoeck 2004), the relationship between th e 3 ice col- 
umn densities can be visualised using a ternary plot dCook et all 
120101) . This is shown in Fig.[2] (d), where column densities are 
shown as a percentage of the total ice column density. Each ver- 
tex of the triangle represents mixtures of pure water ice (magenta), 
80% CO ice (orange) and 80% CO2 (cyan), with the abundance of 
a particular species decreasing with increasing distance to the ver- 
tex. For each YSO, the relative ice abundance of each species can 
be read by following the dashed lines and reading the scale with the 
corresponding colour (see above). 

The LMC and Galactic YSOs exhibit different ice properties 
as they clearly separate in the ternary diagram. As suggested above, 
CO and CO2 vary in a concerted way. Therefore the only way for 
the LMC sample to be made to occupy the same region in this di- 
agram as the Galactic sample is for either the CO and CO2 abun- 
dances to be reduced, or the water column density to be increased 
(roughly by a factor 1.6—2). In other words in the LMC either CO 
and CO2 ices are overabundant or water ice is depleted. In Fig.|2] 
(d) we also show the ice a bundances of low-luminosity G alactic 
YSOs (small green circles. IPontoppidan et ai1l2003l 120081) . Low- 
and high-luminosity (massive) YSOs separate well in this diagram: 
massive YSOs are locate d further away fr om the CO vertex, an 
effect of CO desorption dCooketalj|2010h related to the higher 
temperatures reached in their envelopes. Furthermore, the Galac- 
tic low-luminosity and LMC objects do not occupy the same re- 
gion in the ternary diagram. This cautions ag ainst the similarities i n 
CO2 ice properties found in these samples bv lOliveira et al.l (2009). 
and reinforces the importance of jointly analysing the 3 major ice 
species. To avoid the added complication of evolutionary effects in 
the YSO envelopes, our analysis is based on the comparison of the 
MC sample to massive Galactic YSOs, that are in the sam e evolu- 
tionary stages dOliveira et a l. 2009; S himonishi et al.l 2010). 

None of the SMC sources show CO ice in their spectrum (up- 
per limits are provided in TableQ]and are shown in Fig.|2j. This 
seems to suggest that CO ice does not form or is swiftly desorbed 
around SMC YSOs. However, only for one SMC YSO is the C0 2 
column density above the CO2 threshold (~4.7 xlO 17 molecule 
cm -2 ) for the detection of CO ice (TableQJ. We believe CO ice is 
depleted in the SMC (the upper limits are stringiest in the sources 
with largest water and CO2 column densities) but our data are not 
yet conclusive. In the ternary diagram the SMC sources obviously 
appear at the limit of low CO column densities. 

4 DISCUSSION 

Of the three ice species investigated here, CO is the most 
volatile: desorption temperatures are ~ 20—30, ~ 50—90, and 
~ 100— 120 K, respectively for CO, CO2, and H 2O, depending on 
the ice mixtures in which the species reside dFraser et al.l l200ll ; 
Ivan Broekhuizen et al.ll2006T) . The largest contribution to CO2 ice 
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Figure 3. Schematic of the column densities (not to scale) of the three 
ice species (for a given H2 column) compared between the Galaxy, LMC, 
and SMC. Within the molecular cloud, the outer region (light pink) of the 
YSO envelope contains only water ice (the least volatile species) while the 
shielded interior region (dark pink) contains all three ice species. The main 
observed effects are lower water ice column going from the Galaxy to the 
LMC, and a reduced CO ice column when moving on to the SMC. 



column density arises from chemical reactions within the water- 
rich i c e layer, both in Gala ctic and LMC objects 1 Gerakin es et al.l 
1 19991 ; IOliveiraetalj r2009). On the other hand, the largest con- 
tribution to CO column density comes from the CO-rich layer 
that freezes out on top of t he water-rich lay er, at higher densi- 
ties deeper within the cloud dGibb et ai]|2004l at least for Galac- 
tic YSOs). Thus, CO su rvives only in the more shielded regions 
dHollenbach et alj|2009l) . Consequently, environmental changes in 
the molecular clouds can change the balance of species in the icy 
layers. 

The ternary diagram (Fig.[2](d)) clearly shows that ice abun- 
dances towards YSOs in the LMC are distinct from those towards 
Galactic YSOs. To exp lain the larger N(CQ 2)/ N(H 2 0) ratio ob- 
served in the LM C, bothlOliveira et al.l|2009l) and lShimonishi et all 
d20ld see also lShimonishi e t al. 2008) suggest that CO2 produc- 
tion is enhanced in the LMC. These autho rs propose this is due 
to the stro nger UV field dWeltv et al.l l2006) and higher dust tem- 
peratures dBernard et al.l 120081) observed in the LMC, that have 
been shown in laboratory experiments to increase C O2 produc- 
tion dD'Hendecourt et al.ll986! ; lRuffle & Herbstl200lh . While CO2 
forms in the icy grain mantles, CO forms in the gas-phase and 
freezes out onto cold dust grains. If the temperature of the grain 
surfaces is on average higher as proposed above, CO ice abundance 
would be decreased. Thus it is difficult to see how solid CO abun- 
dance can be enhanced in the harsher ambient conditions in the 
LMC. 

If both CO2 and CO ice abundances appear enhanced towards 
LMC objects but in such a way that relative CO-to-C02 abundances 
are unchanged (see Fig.|2](c)), another possibility is that the wa- 
ter ice column density is reduced towards YSOs in the LMC. In- 
deed, the largest water column density shown in Fig.|2] are ~ 50 
and ~ 100 x 10 17 molecule cm , respectively in the LMC and the 
Galaxy (lower in the SMC), even though this could be a selection 
effect. We explore a scenario (Fig.O to explain the reduction of the 
water ice column density at lower metallicity. 

Water ice forms easily in molecular clouds. However, accu- 
mulation of water ice only occurs once there is enough shield- 
ing and the rate of ice formation overcomes the rate of photo- 
desorption (assuming the dust grains are cold enough and thermal 
desorption is negligible). The Av -threshold for water ice freeze- 
out ( Ayf) is affected by th e strength of the incident radiation 
field dHollenbach e"t~al]|2009l) . varying between Av / ~ 3 mag (for 
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typical Galactic molecular clouds) to Ay / ~ 6 mag. The reduced 
metallicity of the MCs co mpared to the Galax y result s in a lower 
dust:g as ratio (e.g. LMC, lBernardetai1l2008l : SMC, Lerov et al. 
2007); this combined with the harsher UV interstellar radiation 
field (e.g. IWeltv etafl l2006h probably leads to enhanced photo- 
desorption of water and reduced freeze-out in the outer H2 en- 
velopes of molecular clouds in the MCs. Hence the "snowline" 
for water moves inwards in the YSO envelope, reducing the col- 
umn density of water ice as observed in the MCs (Fig.[3j». A re- 
duced availability of water-rich ice is also in agreement with the 
lower fraction of CO 2 locked in the polar component as found by 
IOliveiraetaT] ( l2009h . 

Within the shielded denser regions of the YSO envelopes, CO2 
and CO are abundant in the LMC, and we see no significant dif- 
ference between Galactic and LMC YSOs in terms of these ices. 
CO ice is even more sensitive to shielding than water ice (the Ay- 
threshold for the detection of CO ice in Galactic sources is> 6 mag, 
iBergin et al]|2005l : Iwhittet et alj|2007h . However, gas-phase CO 
also requires high density to accrete onto the icy grains. For the 
same amount of dust (i.e. equivalent s hielding) the correspo nding 
gas density is higher in the LMC (c.f. Ivan Loon et al. 2010b) due 
to the lower dust:gas ratio. These two effects could balance in such 
a way that CO freeze-out still occurs. 

In the Galaxy, the Ay-threshold for the dete ction of CO2 ice 
is ~ 4 mag ( Ber gin et al.l2005l : lwhittet et alj2007l) . consistent with 
that of water ice within the large uncertainties. Our scenario relies 
on an extinction threshold difference between the two species, that 
is not currently detected in the Galaxy. However, it is possible that 
such a difference is more pronounced in the low metallicity/high 
UV environments prevalent in the MCs. Even if CO2 ice is also af- 
fected by the reduced shielding in the outer cold regions, the higher 
density argument put forward in the previous paragraph would also 
favour an enhancement in CO-based CO2 f ormation in the LMC, 
as it is indeed observed dOliveira et alj2009l) . 

In the SMC the CO2 ice abundance seems depressed and no 
CO ice is detected (see caveats discussed in the previous section). 
This could be a direct effect of the lower CO gas-phase abundance 
(e.g. lLerovetaT1l2007l) , and thus a clear metallicity effect. Further- 
more. lvan Loon e t al. (2010b) found that typical dust temperatures 
within YSO envelopes in the SMC are higher than in the LMC. 
Both these effects seem to prevent significant freeze-out of gas- 
phase CO in the SMC. Hence YSOs in the SMC show reduced 
CO2 and CO ice column densities when compared to water ice. 

5 SUMMARY 

We find that the observed column densities of CO2, CO and water 
ice towards LMC YSOs are more consistent with the depletion of 
water ice, not an increased CO2 production as previously proposed. 
This depletion can result from the combined effects of a lower gas- 
to-dust ratio and stronger UV radiation field, that forces the onset 
of water ice freeze-out deeper into the YSO envelope, reducing the 
observed column density. This trend may continue on to the SMC. 

We present the first CO and CO2 ice spectra of massive YSOs 
in the SMC. CO2 ice seems to be depleted and no CO ice is de- 
tected, possibly due to the lower gas-phase CO abundance, and 
higher dust temperatures in the SMC YSOs. 
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